Dense monoenergetic proton beams from chirped laser-plasma interaction 
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Interaction of a frequency-chirped laser pulse with single protons and a hydrogen gas target is 
studied analytically and by means of particle-in-cell simulations, respectively. Feasibility of gener- 
ating ultra-intense (lO'^ particles per bunch) and phase-space coUimated beams of protons (energy 
spread of about 1%) is demonstrated. Phase synchronization of the protons and the laser field, 
guaranteed by the appropriate chirping of the laser pulse, allows the particles to gain sufficient 
kinetic energy (around 250 MeV) required for such applications as hadron cancer therapy, from 
state-of-the-art laser systems of intensities of the order of 10^^ W/cm^. 
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Interaction of high-intensity lasers with solid targets 
has attracted much interest over the past decade, due to 
its potential utilization in laser acceleration of particles 
[ll4l3J . This has given much needed impetus to efforts di- 
rected at replacing conventional accelerators in the near 
future by compact and relatively low-cost devices based 
on an all-optical acceleration mechanism |14| . In particu- 
lar, hadron cancer therapy [15lll6l| may benefit from this 
trend. 

Several regimes are now in existence for laser-plasma 
acceleration [l7|. For laser intensities of 10^^ — 10^^ 
W/cm^ and solid targets with a thickness ranging from 
a few to tens of micrometers, target normal sheath accel- 
eration (TNSA) is the prevailing mechanism. In TNSA, 
a strong quasi-static electric field is induced at the rear 
surface of a thin foil, as a result of emission and accelera- 
tion of the electrons caused by an intense linearly polar- 
ized laser field. Ion acceleration, by the sheath electric 
field, has been extensively studied [ll-llC|. The regime of 
radiation- pressure-dominant acceleration (RPA) , has be- 
come recently accessible by decreasing the target thick- 
ness. Circularly polarized lasers at normal incidence 
have been employed, which suppress electron heating and 
let all particle species co-propagate as a quasi-neutral 
plasma bunch in front of the light wave ("light sail" 
mechanism) [l8| . Despite recent experimental Il9l and 



theoretical [20| progress, clinically useful ion beams [21 1 
have not yet been produced within a scheme which oper- 
ates at currently available laser parameters. 

In this Letter, we demonstrate the theoretical feasi- 
bility of creating proton beams, of unprecedented energy 
and quality, from illuminating a hydrogen gas target with 
an appropriately chirped laser pulse of intensity accessi- 
ble by state-of-the-art laser systems ^]. The basic idea 
of our model stems from the realization that an incom- 
ing highly relativistic laser pulse quickly ionizes hydro- 
gen in the cell and accelerates the electrons away from 
the much heavier protons, as the pulse intensity rises. At 
high enough laser intensities the protons get accelerated 



directly by the laser field. Chirping of the laser pulse en- 
sures optimal phase synchronization among the protons 
and the laser field and leads to efficient proton energy 
gain from such a field. 

We derive an expression for the energy transfer during 
interaction of a particle with chirped unfocused as well 
as pulsed electromagnetic fields. Results obtained ana- 
lytically for the particle's energy gain will be supported 
by further simulations which describe the focused fields 
more accurately than the simple plane-wave model. Our 
two-dimensional (2D) particle-in-cell (PIC) simulations 
reveal that proton beams of energy around 250 MeV, en- 
ergy spread of about 1%, and density of 10^ particles per 
bunch, can be produced. Beams of such quality may po- 
tentially be suitable for use in hadron cancer therapy as- 
suming the laser-plasma interaction can take place close 
to the patient [23[. Following acceleration, the beams 
have to be coUimated and guided 2J] to compensate 
larger variations in the laser pulse and its chirp which 
is still challenging to control at high intensities. Future 
laser systems at the ELI or HiPER facilities j25l, |26[ may 
be utilized to produce monoenergetic multi-GcV proton 
beams. 

A point particle of mass M and charge Q acquires 
relativistic energy and momentum, respectively, oi £ = 
jMc^ and p = '-fMc(3, where /3 is the velocity of the 
particle scaled by c, the speed of light in vacuum, and 
7 = (1 — /3^)~^/^, when interacting with the fields E and 
B of an intense laser pulse. It suffices in many applica- 
tions to represent the fields of the beam by plane waves. 
A plane wave propagating along the z— axis and polarized 
along the x— axis, may be represented by E = xEof and 
B = yEof /c. The dependence in / on the space-time co- 
ordinates is through the combination cot — kz, with cj the 
frequency and k — uj/c the wave number. This leads to 
a break-up of the energy-momentum transfer equations 
(SI units) into four component equations, namely 
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FIG. 1: (color online). The normalized electric field at focus 
of a plane- wave pulse employing Eq. ((5|. The parameters 
used are: A = 1 fim, and r = 50 fs, ^o = and f; = 4s. 



■^hl^z) = aoJo/S^f, — = aoJoPxf, (2) 

where a = QEq/{Mcuiq) is a normalized laser field 
strength and ujq will be defined below. From the sec- 
ond of Eqs. ^ we immediately identify a first constant 
of the motion, namely 7/3j, = ci . Furthermore, the con- 
stant 7(1 — Pz) = 70(1 — Pzi)) — C2 may be arrived at 
by subtracting the second of Eqs. ^ from the first, and 
integrating. Substituting the constants ci and C2 via /3y 
and Pz into 7"^ = I - {^l + P^ + fil), the following key 
expression for the energy of the particle, scaled by its rest 
energy, may be arrived at 
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The combination 77 = a;o(t — z/c) will be used often 
below. Note that with the help of drj/dt — uja{l — (Hz), 
the first of Eqs. (H]) may be formally integrated to give 



iPx = loPxO + a fiv')dr/ 



(4) 



no 



where /S^o = (ix{m) and 70 = (1 - /3g) 2 . 

Frequency chirping amounts to letting lo vary with 
time in some fashion. Experimentally, laser pulses with 
a near-uniform spectral intensity over two octaves are 
feasible [27]. The phase-coherent synthesis of separate 
femtosecond lasers |28| and the recent synthesis of mul- 
tiple optical harmonics [29[ put forward lasers with an 
even wider frequency range. We will work with the lin- 
ear chirp a; = Wo -(- 60 (i ~ z/c) [30-32]. with cjq denoting 
the frequency aX t — Q and z — Q, and 60 having a unit of 
s~^ . The chirped frequency thus becomes to — wo(l+fe^), 
where h = bo/uiQ is a dimensionless chirp parameter. 
Also, dependence of the fields on the space-time coor- 
dinates may be rewritten as tot ~ kz -^ jj + brf' , and the 
chirped field function /(wi — kz) — > f{ri + hrf). We will 
initially work with 



f = cos{(f>Q + ri + br]'^)g{r]); 5(77) = exp 



2s2 



(5) 



where (^o is a constant initial phase, 5(77) a pulse-shape 
function, and s is related to the pulse duration t (full- 
width at half-maximum) via s ~ liJqt 1(2^/2 In 2). In our 
calculations, we choose a shift in 77 in the envelope func- 
tion 3(77) denoted by 77 = 4s. Eq. ^ approximates an 
actual laser pulse in which the non- vanishing pulse inte- 
gral may be compensated by a quasi-static tail. 

To illustrate the mechanism of acceleration, we show in 
Fig. [l]the normalized electric field E/Eq as a function of 
77. Interaction of a particle with the unchirped pulse will 
result in no energy gain, due to the plane- wave symmetry: 
gain from a positive part of the field gets canceled by 
loss to an equally strong negative part. However, the 
chirping breaks this symmetry. Thus, interaction with 
the low- frequency and strong positive part of the chirped 
pulse, which extends over roughly one half of the pulse 
duration, results in net energy gain. 

For a particle initially (77 — rja) at rest, ci = and 
C2 = 7o = 1- Thus, Eqs. ([3|) and ^ give the following 
expression for evolution in 77 of the particle kinetic energy 

i^(77) = (7-l)Mc2 = ^(7/3,)2. (6) 

Using the initial conditions and Eqs. (j4]) and (O, Eq. 
([6|) can be analytically integrated in terms of a complex 
error function. Thus, Eq. (j6|) gives the kinetic energy of 
the particle explicitly, and helps us determine the value 
of the chirp parameter h that maximizes it. Note from 
Eqs. (HI) and ([5]) that the proton energy scales linearly 
with the laser intensity. 

In order to test the applicability of the analytic plane- 
wave model, we compare results based on it with those 
stemming from the use of focused laser fields. For the 
focused fields, the Lax series expressions will be adopted, 
according to which all components E^ , Ey, Ez, By , and 
Bz [33| are given in powers of the diffraction angle 
e = X/ttwo , where wo is the waist radius of the beam at 
focus. Simulations are performed to solve the equations 
of motion numerically for a single proton submitted to 
the laser fields from an initial position of rest at the ori- 
gin of coordinates (770 = 0). In the simulations, wq = 5A 
gives £ = I/Stt <C 1. Hence, terms in the Lax series up to 
0{e^) only need to be retained. The simulations demon- 
strate clearly that the two models predict exit kinetic 
energies that agree to within 2—3% for the parameters 
used in Fig. [2] The linear scaling of the exit kinetic 
energy with the laser intensity given by the plane-wave 
model only holds true for intensities up to about 5 x 10^^ 
W/cm^. For highly relativistic particles the exit kinetic 
energy scales as the square root of the laser intensity, 
which is the case in conventional TNSA [J], too. 

In Fig. [2] (a) and (b) the exit kinetic energy i^ of a 
single proton is displayed as a function of b in the limit 
of large 77 and for optimal 6 as a function of 77/277, respec- 
tively. The energy gain is found to peak globally at the 
chirp parameter values of 6 '-^ —0.003033 and —0.002979 
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FIG. 2: (a) Exit kinetic energy gained by one proton as a 
function of the dimensionless chirp parameter 6 in the hmit 
of large rj. (b) Evolution of the proton kinetic energy with 
rj/2n for optimal b — —0.002979. Peak power of the focused 
laser system is 1 PW, which corresponds to a peak intensity 
of 2.55 X 10^^ W/cm^ for wo = 5A, r = 50 fs and A 
Other parameters used are: (po — 0, and fj — 4s. 
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for the plane-wave and focused pulses, respectively, and 
for the specific laser system parameter set used. Other 
optimal chirp parameter values, which can be determined 
in like fashion, lead to local, less pronounced kinetic en- 
ergy maxima. The fact that the optimal values of b are 
so close (and so are the corresponding exit kinetic ener- 
gies) demonstrates that the analytic solution can be quite 
reliable and that the plane-wave representation closely 
describes the physics involved, at least for the parame- 
ter set used in which wq 3> A. On the other hand, it is 
obvious that dependence upon b of the exit particle ki- 
netic energy is sensitive, especially in the regime where 
high energy is sought. In fact, the plane- wave calcula- 
tion yields the uncertainty SK « l/2{d^K/db'^)\biSb)^, 
which results from an uncertainty Sb ^ 10^^ in deter- 
minig b. The plane-wave calculation yields SK ^ 0.195 
MeV, or about 0.074%. When one takes 6b ^ 10"^ in- 
stead, the uncertainties increase by two orders of magni- 
tude (SK - 19.5 MeV and 7.4%). 

Having studied the single-particle aspects, we now 
present typical 2D3V-PIC (two spatial dimensions and 
three momentum space degrees of freedom) simulation 
results of the interaction of the chirped laser pulse with a 
pre-ionized hydrogen target. The target might be either 
an expanding hydrogen cluster available at sizes rang- 
ing from 1 nm to 1 /J,m 3J] or part of a hydrogen gas 
jet. We are using the following simulation environment: 
The spatial resolution of our simulation box is given by 
Az = Ax = A/100, where the laser wavelength is still 
assumed to be A = 1 fim. The particle number per cell 
is 100, for both protons and electrons. The x-polarized 
laser enters the simulation box from the left and prop- 
agates in the z-direction. For the fields we choose a 
simple modification of the plane wave fields given by 
E = xEofe-^^/^^'""'> and B = y{Eo/c)fe-^^/^^'"»\ with 
/ from Eq. ([5]) and the factor e"^ /(2™o) mimicking spa- 
tial focusing in the polarization direction. The target 
dimensions are assumed to have a length of 0.2 A in the 
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FIG. 3: (color online). Snap-shots (a) of the electron and 
proton density distribution during laser-plasma interaction 
and (b) of the proton density distribution after laser-plasma 
interaction for various times. The laser peak intensity is 2.55 x 
10^^ W/cm^ with Wo = 5A, r = 50 fs and A = 1 nm. 



laser propagation direction, and extension 0.6A in the 
transverse direction. Here, we assume the electrons (in 
the hydrogen gas) to be underdense Ue = 0.1 ric, where 
Tic — 1.1 X 10^^ cm~'^ is the critical density for A = 1/im 
wavelength. For the protons, the total number per shot 
being accelerated as one bunch amounts to « 10^. 

Fig. [3] provides an insight into the acceleration mecha- 
nism and the plasma dynamics. From Fig. [3] (a), we can 
see that when the laser pulse approaches the gas target, 
the electrons get blown off the target. The density pro- 
file of the electrons follows approximately the laser field 
oscillation. At this time (t — 16 laser cycles) the proton 
distribution almost maintains its initial shape. After the 
interaction of the gas target with the laser pulse, we see 
from Fig. [3] (b) that the protons are accelerated as a 
symmetric bunch of homogeneous density (average den- 
sity at i = 50 laser cycles: rip « 0.04 ric). Note that the 
proton bunch is emitted at 70° relative to the laser prop- 
agation direction (z-direction). Due to the comparatively 
low density of the plasma target and, hence, suppressed 
Coulomb explosion, the time-dependent beam divergence 
is low. For a macroscopic distance of 30 cm, the cross- 
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FIG. 4: (color online). Proton kinetic energy distributions 
after the interaction of a chirped laser pulse (6 — —0.003033) 
with a hydrogen gas target (left peak) and with an ensemble 
of protons without particle-particle interaction effects (right 
peak). The laser peak intensity is 2.55 x 10^^ W/cm^ for both 



pulse with an underdense hydrogen gas target. The re- 
quired laser peak intensity of about 10^^ W/cm^ is within 
the range of state-of-the-art high-intensity laser systems 
[22| . Furthermore, we developed an analytical model for 
determining the optimal pulse shape and the final kinetic 
energy of the particles, which is in good agreement with 
the performed 2D-PIC simulations. Our work has mainly 
been concerned with proton acceleration to energies and 
densities required for hadron therapy. Following accel- 
eration, ion beam shaping 2J] beyond the scope of the 
present Letter has to be applied. 
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is supported by the Arab Fund for Economic and Social 
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Alhance HA216/EMMI. 



sectional area of the beam is about 0.98 cm^. This num- 
ber roughly corresponds to a divergence angle (opening 
angle of the cone defined by the beam) of only 2°. From 
the particle dynamics it is obvious that the prevailing 
part of the proton's kinetic energy is transferred via di- 
rect interaction with the laser field. 

In Fig. 2] we compare the exit kinetic energy dis- 
tribution of 3000 initially randomly distributed non- 
interacting particles (cf. [35|, |36| for the method) in a 
spatial volume with the same dimensions as the gas tar- 
get (cf. Fig 13] (a)), with the one stemming from the laser 
- gas target interaction. In the case of pure vacuum 
acceleration K = 258.3 MeV ± 1.2%, while for the laser- 
plasma-cell acceleration K = 245.2 MeV ± 0.8%. The 
discrepancy between the reported mean kinetic energies 
is about 5%, which can be attributed to particle-particle 
interaction effects. Moreover, the simulation of the ran- 
domly distributed non-interacting particle ensemble has 
been carried out in three spatial dimensions, whereas the 
PIC simulation is two dimensional. Enlarging the initial 
plasma distribution leads to an increase in the energy 
spread and to a decrease in the mean kinetic energy of 
the created proton beam. For example a gas target with 
initial length of lA (5 times larger than before) yields 
K ~ 245.1 MeV ± 5.4%. However, employing some ve- 
locity filter one could maintain the beam quality. In order 
to reach similar particle energies by using an unchirped 
laser pulse one would have to raise the laser intensity by 
three orders of magnitude, to lO^^'W/cm^. At such an 
intensity the chirped laser scheme (employing the same 
parameters as in Fig. [2]) already yields monoenergetic 
proton beams with K = 17.3 GeV ± 1.0%. 

In summary, we have demonstrated the theoretical fea- 
sibility of creating a dense proton beam (10^ protons per 
bunch) of high energy (w 250 MeV) and good quality (en- 
ergy spread ^ 1%), from the interaction of a chirped laser 
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